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Abstract 



It seems widely accepted by the turbulence community that the inter- 
mittency observed in fully turbulent flows is closely related to the exis- 
tence of intense vorticity events, localized in time and space, also known 
as coherent structures. We describe here an experimental technique based 
on the acoustic scattering phenomenon allowing the direct probing of the 
vorticity field in a turbulent flow. In addition, as in any scattering exper- 
iment, the information is in the Fourier domain : the scattered pressure 
signal is a direct image of the time evolution of a well specified spatial 
Fourier mode of the vorticity field. Using time-frequency distributions, 
recently introduced in signal analysis theory, for the analysis of the scat- 
tered acoustic signals, we show how the legibility of these signals is signif- 
icantly improved (time resolved spectroscopy). The method is illustrated 
on data extracted from a highly turbulent jet flow : discrete vorticity 
events are clearly evidenced. The definition of a generalized time-scale cor- 
relation function allows the measurement of the spatial correlation length 
of these events and reveals a time continuous transfer from large scale 
towards smaller scales (turbulent cascade) . We claim that the recourse to 
time-frequency distributions leads to an operational definition of coherent 
structures associated with phase stationarity in the time-frequency plane. 
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Intermittency |IJ is usually revealed by analyzing the statistical properties 
of the flow deduced from single point measurements (e.g. the time record of the 
flow velocity at one point obtained by hot wire anemometry). In the past few 
years, the turbulence community has reached a consensus about an acceptable 
definition of intermittency, this latter being defined as deviations (or corrections) 
from the power law scaling of the longitudinal velocity structure functions in 
the inertial range: S p (l) — ((5v\\(l)) p ) oc £«> initially proposed by Kolmogorov 
in 1941 (K41 theory). Indeed, the K41 theory predicts that Cj> = §> whereas 
clear deviations from this linear law are experimentally observed Recent 
results suggest that intermittency is closely related to the presence of localized 
(in time and space) vortices in the turbulent flow. Such objects, presumably 
with a filamentary shape, have been observed in some experiments |J and 
numerical simulations J^, [|. Theoretical models [Q assuming a hierarchy of 
coherent structures in the form of vortex filaments, predict a relation for the 
C p = f(p) in very good agreement with experimental data. However, a clear 
definition of coherent vorticity structures is still lacking. 



There exists a vast amount of theoretical work dealing with the acoustic 
scattering phenomenon by velocity fields. Most of it was initiated nearly fifty 
years ago and started with the papers authored by A.M. Obukhov R.H. 
Kraichnan 0, Chu and Kovasznay p0| , to mention but a few. In a more re- 
cent paper, F. Lund |]ll"|| has established, under reasonable assumptions, a linear 
relation between the scattered amplitude of a plane acoustic wave incident on 
a turbulent flow and the spatial Fourier transform of the vorticity field. The 
purpose of this paper is to demonstrate how one can take advantage of the spec- 
tral nature of the measurement process, in an acoustic scattering experiment, 
for the detection of spatially and temporally localized vorticity events. Indeed, 
there are numerous examples in physics of disordered media, where scattering 
techniques (e.g. light and neutron scattering in condensed matter physics) pro- 
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vide useful and pertinent informations about the statistical properties of these 
systems. 

This experimental work aims at introducing, on experimental basis, a new 
method for investigating turbulent flows, serving as an alternate to usual single 
point measurements like hot wire anemometry. We place the emphasis on some 
experimental evidence of the existence of coherent and localized (i.e. intermit- 
tent in both space and time) structures. Using a dual channel setup (acoustic 
interferometer), we are able to experimentally demonstrate the effectiveness of 
the length scale selection (through the selection of a well defined spatial scatter- 
ing wave vector) involved in acoustic scattering measurements. Further, as in- 
tcrmittency implies non stationarity of physical quantities (think about a vortex 
as being localized in both time and frequency domains, appearing or disappear- 
ing within the analyzed volume), the second part of this paper focuses on non 
stationary spectral analysis. The proposed approach is based on the recently in- 
troduced |L2j Reduced Interference time-frequency Distribution (RID), of which 
the properties are briefly discussed. RID analysis is then applied to signals from 
an experiment of acoustic scattering by a fully turbulent jet flow. On a single 
RID reprentation of the scattered pressure signals, one clearly oberves well lo- 
calized energy packets at a given spatial wavevector. In order to get further 
insight on these objets, we compare the RID representations at two different 
wavevectors simultaneously recorded. This dual representation demonstrates a 
significant correlation of the time occurences and of the velocities (correlated 
phase derivatives due to correlated doppler shifts). As this time- frequency cor- 
relation decreases with the separation between the two wavevectors, we conclude 
that the different patterns selected at different wavectors belong to what can be 
defined as a coherent structure. Notice here that our approach closely matches 
the definition of intermittency formerly proposed by G.K. Batchelor in 1953 
[ p"3| : "The inference, then, is that there is an uneven distribution, in space of 
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the energy associated with the large wave-number components of the turbulence, 
and that the higher the wave-number, the more does the associated energy tend 
to occur in confined regions of space ( meaning that if a Fourier resolution of the 
velocity field within a region of activity were made, the amplitude of the com- 
ponent at the relevant wave-number would be found to be large, while a Fourier 
resolution of the velocity field within a region of quiescence would give a very 
small amplitude; the amplitude of the component for the field as a whole will 
lie somewhere between these two amplitudes)" . The paper is organized as fol- 
lows : part 1 is devoted to the presentation of some theoretical aspects and 
assumptions underlying scattering of acoustic waves by a turbulent flow field 
with a particular emphasis on the role played by the vorticity field in this inter- 
action. In part 2, we describe the experimental setup as well as the statistical 
properties of the turbulent jet flow under investigation. The data acquisition 
procedure and the preprocessing of the recorded data are detailed in part 3. 
Part 4 deals with some theoretical issues of time frequency analysis and the role 
of the ambiguity function for introducing the RID. Finally, experimental results 
are displayed and discussed in part 5. 

1 Acoustic scattering by a turbulent flow 

As quoted by Nelkin fl4| , vorticity is intrinsically a non-local quantity, and 
therefore it can hardly be measured by using local probes (e.g. hot wire anemome- 
ters) |p^| . Although vorticity could " conceptually" be measured by local probes, 
the main problem resides in that one needs at least two probes, separated by 
a finite distance in order to compute the gradients involved by the curl opera- 
tor. From an experimental point of view the situation is even worse: at least 
nine probes are needed to reconstruct the three components of the vorticity, 
within the appropriate Taylor Hypothesis (see p6[). Since these probes have a 
finite extension and finite time response, the computation of the gradients at 
the smallest scale of the turbulent flow becomes problematic |i~7|, [l8| . It is worth 
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noticing here that an alternate experimental has been proposed by Noullez & 
al. |lE{ ] in order to probe directly the spatial transverse velocity increments in 
a turbulent flow. 



Taking these latter remarks into account, one can find several theoretical 
[§, §, fo|, experimental @ §l[ ||, §| §|[ |f| and numerical works (if] address- 
ing the direct measurement of the vorticity field using acoustic scattering. As 
in any scattering experiment, one expects to be sensitive to the spatial Fourier 
distribution of the scatterers (as the amplitude of the scattered wave is ex- 
pressed in terms of some coherent average on "secondary sources" ) over a finite 
volume defined as being the cross section of the acoustic beams (incident and 
scattered). In a recent paper, Lund et al. ]Tl|] derive a linear relation between 
the Fourier transform of the scattered ultrasound amplitude and the space-time 
Fourier transform of the vorticity field. In analogy with the more usual light 
scattering phenomenon, the physical mechanism at the origin of acoustic scat- 
tering by vorticity can be thought of as follows : an acoustic wave impinging 
on a vorticity distribution induces fluctuations of the vorticity at the incoming 
sound frequency (by virtue of the Kelvin circulation theorem). Each scatterer 
(vortex) acting as a secondary source will, in turn, radiate a sound wave. The 
coherent average (taking into account the relative positions of the individual 
vorticity elements) over the scatterers distribution results in the emission, out- 
side the vorticity domain, of a scattered acoustic wave. It is important to 
emphasize here that many vortices interact simultaneously with the wave front. 
Such interactions have been numerically investigated in Jl9| in the case where 
a unique vortex is present. In a scattering process, one is sensitive to spatial 
details of the scatterers distribution at a typical length scale measured by the 
wavelength. Note that whereas a light scattering process is usually linear, the 
acoustic scattering phenomenon depicted here stems from the non-linear term 
of the Navier-Stokes equation, and requires much explanations and details com- 
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putations that are beyond the scope of this paper. In particular, one can find in 
[ p"l| a clear explanation of the respective contributions to the acoustic scattering 
cross-section of the vorticity field on one hand, and of the irrotational velocity 
field, induced by the vorticity, on the other. 

Using a Born approximation, Lund et al. obtain the following linear rela- 
tion between the scattered acoustic pressure amplitude and the spatial Fourier 
transform of the vorticity field: 

VscaM 2 . -CQS(6 S ) ve™ D / c ^ A ft 

= 71 l i-co S (e s ) *d {n ^ n(ff ~* " ~ Vo) (1) 

where A and "." stand for the vector product and the scalar product respectively, 
and 

27r / -s 47W o . ,0 S r-n 

qscat = — [yr - v n) ~ sin - - — — lor v ~ v a (2) 

c c 2 \r — 7i | 

In these formula, 

• Qscat is the scattering wave vector 

• n and f are the impinging and scattered wave directions, respectively 

• v a is the frequency of the incoming sound 

• v is the frequency of the scattered sound 

• S is the scattering angle 

• Pine stands for the pressure amplitude of the incoming sound wave (as- 
sumed to be plane and monochromatic) 

• Pscat is the amplitude of the scattered sound wave 

• D and c stand for the acoustical path between the measurement area and 
the detector, and the sound velocity respectively. 

• ft and r are the unit vectors in the direction of the probing and detected 
sound waves 
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It is important to emphasize here that among the hypothesis used in |Tl|] the 
characteristic time of any process within the flow has to be very large compared 
with the period of the analyzing sound wave. This latter remark implies that 
v ~ v Q and that this method is relevant only for flows with low Mach number. 
Equation ([!]) shows that, using acoustic scattering, it is possible to probe one 
component (perpendicular to the scattering plane defined by the vectors n and 
r) at a chosen length scale through the selection of a known scattering wave 
vector q sca t- 

In a previous similar experiment p4[ (see also [pq] ) performed on the regular 
vortex street behind a cylinder at low Reynolds numbers, we have demonstrated 
that the scattering cross section is indeed strongly dependent on the direction of 
the vorticity field accordingly to Equation (1) of the present paper. This latter 
experiment evidences the existence of direct interactions between vorticity and 
acoustic waves : among the nine components of the stress tensor only one vector 
can be defined without ambiguity, i.e. the vorticity (antisymmetric part of the 
tensor). Furthermore, a nice feature of this method lies in the fact that it is a 
non perturbating technique. 

2 Experimental setup 

The acoustic waves are produced and detected by Sell-type transducers (see 
e.g. [^9|), the dimensions of which are 15cm x 15cm. The advantage of dealing 
with such large transducers is to minimize diffraction effects and therefore to 
enable a sharp selectivity of the scattering wave vectors (q sca t defined above), 
both in modulus and direction pi] . As shown on figure (Q), the direction of 
Qscat is aligned with the mean-flow velocity in order to maximize the Doppler 
shift. The magnitude of q sca t, and thus the probed length-scale in the flow, is 
then selected by tuning the incoming sound frequency v Ql the scattering angle 
being held fixed. Indeed, for practical purposes, it is much simpler to hold the 
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geometry of the experimental setup constant. Therefore, 9 S is constant and 
Qscat has constant direction. The scattering angle 9 S was set to 60 degrees, as 
a trade-off between optimizing the sensitivity of the measurements and limit- 
ing the spurious effects induced by diffraction side-lobes pi) , ft can then be 
seen from (|f|), that the Fourier transform (FT) of the scattered pressure signal 
is proportional to the band-pass spatial frequency filtered vorticity field. The 
central frequency of the latter filter is q SC at{vo) = t|<fscat(^o)||- Its width (Aq) 
can be obtained from the diffraction pattern of the Sell transducers [|4j. A 
rough estimate of Aq is given by differentiating q sca t with respect to 9 S , which 
leads to Aq = cos(9 s /2).2A9, where A9 is given by a diffraction ansatz: 
A9 ~ ^j- , L being the characteristic size of the transducer . The volume 
which is analyzed is defined by the cross section of the incident beam and the 
main lobe of the receiver (antenna beam) . Its dimensions in our experimental 
setup were of the same order of magnitude (a bit larger) than the transducers 
dimensions. 

The flow which we investigate in this paper, consists in an axi-symmetric tur- 
bulent jet flow (R\ ~ 600). emerging from a circular nozzle of diameter 5 cm, 
at a high velocity (about 60 ms^ 1 ). Details about our wind tunnel facility can 
be found elsewhere § . The volume of interest was set at a distance of about 50 
nozzle diameters downstream ensuring that the turbulence is fully developed. 
A study based on classical ancmomctry (performed in the same area of interest) 
leads to the estimated following parameters : the average longitudinal velocity 
is about 5.3 ms _1 and the turbulence ratio is about 28% (^6). The third order 
structure function exhibits a clear power law scaling over 1.5 decades (from the 
Taylor scale A = 6 mm up to the integral length scale l = .25 m), with an ex- 
ponent close to f in accordance to the Kolmogorov four-fifths law JlJ . From this 
inertial range scaling, we deduce the value of the mean energy dissipation per 
unit mass e = 13 m 2 s -3 , the Kolmogorov micro scale rj = 131 /xm and the large 
eddy turnover time t a — 1 68 ms. The analyzing sound frequency can be varied 
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from 5kHz up to 200kHz (where sound attenuation effects cannot be neglected). 
Thus we are able to probe a large part of the inertial range, as well as length 
scales belonging to the beginning of the dissipation range p6[ . However in this 
paper, we shall restrict our analysis to a single length scale pertaining to the 
inertial range ( g y ' = 0.6 cm -1 corresponding to a typical length scale of 1.7 
cm). 

3 Data acquisition and pre-processing 
3.1 Single channel measurements 

An experiment consists in measuring and recording time series of the acous- 
tic pressure signal received on a Sell- type transducer (see above). The emitted 
sound wave can be considered as being a monochromatic plane wave (if diffrac- 
tion effects are neglected at first sight). The sampling device (HP3565) has 
16 bit resolution and its sampling frequency is 262144 Hz. Each recording 
channel is provided with its own local oscillator (tuned at v a ), thus enabling 
a numerical heterodyne detection of the scattered signal around the frequency 
of the impinging sound wave. Note also that, in order to reduce phase noise 
along the whole measurement line, the waveform generators (sound emission), 
the sampling clock and the local oscillators (demodulation) are locked on the 
same 10 MHz master clock. Two time series associated with synchronous phase 
and quadrature signals are recorded at a rate of 8192 samples per second, from 
which the analytical complex signal is computed (characterized by a phase and 
an amplitude). Indeed, thanks to the linearity of our ultrasonic detectors, the 
phase information of the scattered pressure wave is preserved (the phase refer- 
ence being that of the incoming sound wave). 

As indicated by equation ([!]) (expressing a convolution in the time Fourier do- 
main), the scattered pressure signal is obtained as a time modulation of the 
incoming pressure signal (carrier wave) by the low frequency signal associated 
with the vorticity of the flow under investigation (for small enough Mach num- 
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bers). One thus expects a narrow band scattered signal with frequencies close 
to the incoming sound frequency v (see figure |^ a). Then, the purpose of the 
above mentioned heterodyne detection is to significantly reduce the amount of 
stored data. Equations (|l|,^) lead (for v ~ v a ) to express the Fourier transform 
of the detected pressure signal as 

cos(fl a )coB(^) ^ 

Vscat\V) CX il±{q S cat, V - V ) (3) 

Qscat 

where q SC at and Qj_ stand for the magnitude of the scattering wave- vector and 
the vorticity component orthogonal to the scattering plane, respectively. We 
have represented on Figure (a) a typical averaged power spectral density of 
the scattered pressure signal (after demodulation) obtained in the turbulent jet 
flow. Notice that in this figure and throughout this paper, the frequency axis 
represents the frequency shift with respect to the incoming sound frequency v Q . 
The spectrum is clearly non-symmetric with respect to the null frequency. This 
latter fact is a simple consequence of the advection of the vorticity field by the 
flow velocity and to our decision to align the scattering wave vector with the 
direction of the mean flow velocity V mean . Indeed, the scattering of waves by 
a moving structure gives rise to a frequency shift Av as a consequence of a 
Doppler effect || : 

Av = V - V Q = —qscat -Vmean (4) 

Notice that this latter relation can be obtained by applying a Galilean trans- 
form in the real space followed by a double Fourier transform on the space and 
time variables. 

From (||) one expects a divergence of the scattered signal around 9 S = 0, for 
which q sca t — and thus Av — 0. Indeed, the existence of direct paths between 
the emitter and receiver (associated to diffraction effects due to the finite size of 
the transducers) results in the presence of a very strong peak in the spectrum 
p scat {Av) around Av = 0. From (^), we are able to relate q sca t to Ais, and thus 
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get the following expression for the received acoustic signal 

Pscat(vo + Al/) OC -^ti±(q SC at, Av) (5) 

Bearing in mind that the scattered pressure signal is recorded after a numerical 
heterodyning procedure that shifts its analytical frequency spectrum around the 
null frequency (the frequency of local oscillator in the lock-in is set to v ), one 
gets 

P S eat rec {Av) OC -^—&±(q S caU Al/) (6) 

Av 

where P sca t T . ec stand for the scattered pressure recorded signal. Taking the 
inverse FT of equation (^J) leads to 

/+oo ^ 
n±(q,cat,t)dt (7) 
-oc 

The preceding equation expresses the recorded scattered signal (after hetero- 
dyning) as being the integral of fij_ (q S cat , t) with respect to the time variable 
As a consequence, the recorded data must be differentiated to give a correct 
representation of Q±(q SC at, t). It is important to notice here that the latter 
differentiation must be performed at a low frequency, since the heterodyne de- 
tection has already been applied. The averaged power density spectrum (PDS) 
of p S cat(t) recorded for the jet flow (see preceding section) and v = 32kHz, is 
shown on figure (|| a) . Figure (|2| b) illustrates the efficiency of the difference filter 
applied to the time series for extracting a PDS proportional to d,±(q scat , Av). It 
is worth noticing at this point that this differentiation ansatz is nothing else than 
a transposition of the Taylor hypothesis (heavily used in single point measure- 
ment). Actually, it suffers from the same limitations with regards to the level 
of turbulence. However, as explained below, these limitations can be overcome 
by the resort to the time-frequency distributions. 

3.2 Two-channels measurements 

According to (|2|), a given spatial spectral component of the vorticity field may 
be probed by using different sets (9 s ,vq), associated with a couple of trans- 
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ducers (one emitter, one receiver) which will be referred to as a channel in 
the rest of the paper. The acoustical paths corresponding to the null scatter- 
ing angle for which divergence occurs (thus explaining the presence of a strong 
peak at zero frequency shift, see ([|)) is then different for each set of ultrasonic 
emitter-receiver. As a consequence, it is expected that the parts of the signals 
corresponding to the null scattering angle, recorded on each channel are mutu- 
ally non coherent. Following the same reasoning, it is expected to find a high 
degree of coherence between parts of signals that are due to scattering effects 
from the same measurement area (cf preceding sections), provided that both 
channels are tuned to probe the same scattering wave-vector: 

4^,1 . ,08,is 4tw ,2 . ,08,2s , Q s 

sm( — ) = - sin(^-) (8) 

c 2 c 2 

where fo t i,0s,i stands for the ultrasound frequency and scattering angle associ- 
ated with channel i. 

The additional hypothesis in the 2-channcl measurement scheme, is that the 
measurement areas defined by the 2 sets of transducers are actually identical. 
This point is crucial in our experiment, as any mismatch in the measurement 
areas would preclude any possibility to compare scattering signals, which would 
then be associated to different spatial distributions. The setting of the exper- 
iment is derived by carefully aligning the transducers with laser diode rays. 
Then, the lengths of the acoustical paths for each channels are set to the same 
value (up to a precision of lmm)Q in order to insure equivalent propagations 
delays on both channels. This is achieved by checking the phase of the cross 
channels inter-spectrum estimates, and shifting the transducers along its main 
diffraction axis accordingly: for null propagation delay, the slope of the phase 
with respect to frequency must be zero, as group delays are expected to be 
identical on both channels. 

Figure (||) illustrates the fact that the divergences at zero-frequency shifts, mea- 

1 The propagation time of acoustic waves over this distance is about 3.10 - e s, i.e. much 
smaller that the sampling period. 
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sured on two channels verifying (ra) are mutually non coherent, whereas the part 
of the spectrum corresponding to the scattered waves from the same measure- 
ment area (Doppler shifted) are highly similar (high cross-coherence values). 
Furthermore, the values obtained for the cross-coherence do not depend on 
whether the differentiation (see preceding paragraph) was performed, as the 
same linear filter (time differentiation) is applied to both channels. This latter 
observation illustrates the pertinence of the differentiation for extracting the 
interesting part (i.e. the signature of the vorticity) in the recorded signal. The 
phase of the cross-channels inter-spectrum is shown on the last plot of figure 
: one sees that the phase slope takes low values within the frequency domain 
where the spectrum of scattered signals is significant. 

So far, we have dealt with second order statistical properties of the turbulent 
flow, and considered long enough records, to insure that the stationary hy- 
pothesis was correctly met. However, the spectrum of the recorded time-series 
corresponding to one spatial frequency component of the vorticity field, turns 
out to have an important width (Equivalent lOdB bandwidth= 400Hz). This 
may be ascribed to either fluctuations of the advection velocity (thus leading to 
fluctuations of the Doppler shift), or to finite duration of vorticity structures. 
In this latter case, one then expects to record "sine-wave trains' of finite du- 
ration. Notice that both effects are likely to occur simultaneously. Therefore, 
one is naturally led to look for an analyzing tool allowing a joint representation 
in time (in order to access "vortex life time") and frequency (for evidencing 
Doppler shifts, or cquivalently, advection velocity fluctuations). 

4 Time-frequency analysis 

When finite duration events of vorticity are present, the scattered acoustic signal 
will exhibit a time- varying spectrum, as the frequency of the measured pressure 
signal drops from vq to vq + 8v when a vorticity event exists, that has energy 
in the spatial Fourier domain at the wave- vector q sca t- Due to the randomness 
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of the velocity field which advects the vorticity field, the finer details of the 
time dynamics of individual events are wiped out when average spectra of the 
scattered signals are computed. Indeed, on the average, the scattered spectrum 
has a Gaussian shape reflecting the statistical properties of the flow velocity 
[ ^6| . Notice that an infinitely thin and localized vortex (theoretical) would have 
contributions at all wave-vectors in the spatial Fourier domain. 

Our aim is thus twofold : first, detecting presence/ absence of such events in 
the recorded signal, and second, estimating the frequency shift Av at which this 
occurs, in order to characterize its relative velocity. We are thus faced with the 
need for a representation of the signal which preserve simultaneously the time 
and frequency information. 

Temporal fluctuations of the spectral distribution can be approached by the 
common spectrogram analysis. Watkins among others, has addressed this 
problem and provides examples of reasonable and unreasonable spectrogram 
results in the analysis of bio-acoustics signals (marine mammal sounds). How- 
ever, the spectrogram involves a sliding time window which intends to capture a 
portion of the signal which is sufficiently restricted in time so that stationarity 
assumptions are approximately met. Furthermore, the presence of the window 
results in a time-frequency distribution (TFD) exhibiting both temporal and 
spectral leakage. 

The Wigner-Ville distribution (WVD) avoids the problems of windowing the 
signal but the presence of interference terms between signal components, due to 
the quadratic nature of the WVD (see next paragraph) |38| often precludes 
its applicability. Recent studies have shown that a class of smoothed WVD, 
the Cohen's Class of TFDs, allows to reduce the amplitude of the interference 
terms, while preserving covariance properties of the WVD in both the time and 
frequency domains. Cohen |35| and Flandrin J3?j have recently provided 
excellent reviews of TF analysis, including recent developments. It is worth 
noticing here that the above mentioned signal representations do not rule out 
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the Gabor-Heinscnbcrg inequalities (see p6fl). A brief presentation of the main 
issues on Cohen's class TF analysis and the key role played by the ambiguity 
function is given below. 

In this paper, we focus on one distribution from the Cohen's class, the Reduced- 
Interference Distribution (RID). It was first introduced by Choi and Williams 
|33), and has the additional feature of preserving the time and frequency sup- 
port of signal components : {RID{t, /) is non zero if the signal is non zero at 
t and has a spectral component at /). These are highly desirable properties if 
one wishes to infer statistical properties of the vorticity distributions from this 
TF analysis. The examples shown were computed with a binomial RID, which 
is a good discrete approximation of the continuous distribution ]39| . 

4.1 The WVD and Ambiguity function 

The Wigner-Ville distribution of z(t) is defined as the Fourier transform (F) of 

R z (t,r) = z (t+ T -)z*{t- T -) 

with respect to the lag variable r. 

W z {t,v) = Fr [z(t+^)z*(t- I)] = F T [R z (t,r)] (9) 

where z(t) is the time signal, z* is its complex conjugate. The operator Ft \ 
denotes the Fourier transform operator, with respect to the variable (.); this 
notation will be adopted in the rest of that paper. From an experimental and 
heuristic point of view, the effect of the operator R z (t, r) is equivalent to a local 
phase conjugation of the signal which enhances the phase derivative related 
to the Doppler shift (P. Flandrin : private communication). Similarly, but a 
with different physical meaning, the symmetrical^ ambiguity function (AF) is 

2 Ordinarily, the AF is defined as the FT of R z (t, r), with respect to t. The symmetric AF 
used here is defined as the inverse FT of R ^lt. r). with respect to t. This does not induce any 
change in the interpretation of the AF (see |35[ ) , and is only a matter of convention. However, 
it allows to r elate the AF to the WVD via a two-dimensional FT with respect to t and r (see 
equation hlh. 
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defined as the inverse Fourier transform (F x ) of R z (t,r) with respect to the 
first variable^] 

A z (6,r) = FT 1 [z(t + l>*(t \)\ = FT l [R*(t,T)] (10) 

Thus, W z (t,v) and A z (9, r) are related by the two-dimensional Fourier trans- 
form 

W z (t,v)= [ [ A z {e,T)e-^ te+VT) d0dT (11) 



4.2 Cohen's Class of Distributions 

Cohen's Class of distributions is defined as 

C z {t,v,<j>) = j J JetM-W-^(6,T)z(t + ^)z'(t-^)dZdrdB (12) 

where <f)(6, r) is the kernel of the distribution^]. 

These relationship can be combined with Eq. @ to show that C z (t,v,4>) 
may be expressed as: 

C,{t,v,4>)= [ [ cf>(6,T)A z (6,T)e-^ te +^d8dT (13) 



Thus, while W z (t, v) is obtained from the symmetric ambiguity function by 
means of two-dimensional Fourier transform, any member of Cohen's Class of 
distributions may be found by first multiplying its kernel 4>{6,t), by the sym- 
metric ambiguity function and then carrying out the two-dimensional Fourier 
transform. The spectrogram and the WVD can easily be shown to be members 
of Cohen's class by choice of the proper kernel in each case. 
The generalized ambiguity function (()(9,t)A z (9,t) is a key-concept in t-f anal- 
ysis which aids one in clearly seeing the effect of the kernel in determining 
C z (t, v, <p). In the ambiguity domain, application of the kernel of Cohen's class 
involves multiplication of the symmetric ambiguity function by the kernel. This 

3 The ambiguity function reduces to a deterministic correlation function in time if 8 is set 
to zero. Similarly, as can be seen in a dual form which starts with Z(u), it can be seen as a 
deterministic correlation of spectra if r is set to zero. 

4 The range of integrals is from — oo to oo throughout this paper. 
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is roughly analogous to application of filters by multiplication of the Fourier 
transform of the signal by the transfer function in the frequency domain. This 
operation manifests itself as a convolution in the time domain, thus, 

C z (t,p,<f>) = F T [ip(t,r) * t R(t,r)] (14) 

where ip(t,r) = Fg [(/)(9,t)]. Similarly, application of the kernel in the TF 
domain involves convolution in time and frequency in Cohen's class of distribu- 
tions. Conventional filtering concepts can be adapted to Reduced Interference 
Distribution (RID) design. 

4.3 The RID approach 

New distributions with reduced interference (RIDs) J3^, |l2[ ||{| have been de- 
veloped in recent years, serving as attractive alternates for the spectrogram and 
the WVD. The RID represents an attempt to preserve most of the desirable fea- 
tures of the WVD while alleviating most of the confusing interference between 
signal components. 

RID kernels exhibit low pass characteristics in both the time-shift (r) and 
frequency shift (9) dimensions in the ambiguity plane, away from the axes. It 
has been shown |l^, [59| that any real valued symmetric, unit area function 
h(t) which is limited to insure time and frequency support properties, such that 
h(t) = for \t\ > h and which has a H{v) \ u —q = 1 is a RID kernel. One simply 
replaces H (y) with H(9t) in order to have a RID kernel in the ambiguity domain 
form. Usually interference terms are located far from the 9 — r origin, so the 
filter can often be designed to provide strong attenuation of these terms, while 
preserving the desirable auto-terms of the individual signal components. The 
form of the RID kernel for the (t,r) domain application is ip(t,r) = j^h(^) 
and computation is usually most convenient in this domain. Preservation of the 
time and frequency marginals is one of the important properties of RIDs. One 
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can recover the energy spectrum (by integrating with respect to time) or the 
instantaneous power (by integrating with respect to frequency). This is not the 
case for spectrograms and many other popular TF distributions (see e.g. p0|). 

5 Application of TF distribution analysis to tur- 
bulence signals 

A nice discrete kernel which has a binomial form has been used in our applica- 
tions. This kernel is computationally efficient and also exhibits desirable char- 
acteristics in terms of TF distributions |^2| . The time frequency distribution 
obtained by applying the binomial kernel takes the following discrete form: 

TF ridM (n, v) = (15) 

r=oo w=+|t| ( \ / i I \ 

T = -OC V = -\T\ V 7 

where z[n\,n — 1 . . . N depicts the sampled time-series under analysis. Func- 
tions h(r) and g(v) are frequency smoothing window and time smoothing win- 
dow respectively. h(0) and G(0) = TF v [g(v)] are forced to one. Note that these 
windows were set to h(r) = rect(N) and g{y) — S(v) in the original paper by 
Williams and Jeong, where rect(TV) stands for the TV-points rectangular win- 
dow, and 5 stands for the discrete Dirac distribution. The binomial TFD is a 
good discrete time approximation of the continuous RID. The binomial TFD 
(hereafter called the tf-rid) was computed for acoustic signals obtained in the 
turbulence experiment depicted in the previous sections. 

In the cases which are presented, two channels were used to probe the same 
measurement area : the experimental setup is carefully arranged in such a way 
that the measurement areas defined by the intersection of the acoustic beams of 
each channel are actually superimposed. This latter setting will be maintained 
for all the experiments presented in the remainder of this paper. First and sec- 
ond channel had the following characteristics (0 Si i = 60 degrees, ^o,i = 20kHz) 
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and (9 S ,2 = 40 degrees, 1/0,2)- Applying eq. (0) and with the transducers ar- 
ranged as described on figure ([!]), one shows that the probed scattering wave 
vectors are co-axial with the mean flow, and oriented downstream. These latter 
channel characteristics correspond to |<? sca t| = 369.6m -1 i.e. the corresponding 
scale (wavelength) is 1.7 cm, well within the inertial range of the flow. The same 
wave vectors are probed onto two different channels by tuning ^ .2 such that 
the condition expressed in eq.(|^) is verified. This latter condition is met here 
for ^0.2 = 29.25kHz. Figure shows the tf-rid obtained for the case where 
the preceding tuning condition is approximately met (actually, 1/0,2 = 29400 Hz 
for this experiment). Smoothing windows h and g were both Hanning windows 
of respective length Lh — 129 and L g = 3. 

The frequency marginals, obtained by integrating the tf-rid over time are shown 
on the left hand side of the plot. This frequency marginal matches the power 
density spectrum (PDS) (an approximation is introduced by the use of smooth- 
ing windows h and g). The center of mass Av cm of the PSD is indicated. 
Replacing Av by Ai/ cm in eq.(Q) allows one to estimate the average (estimated 
over 0.6 seconds) axial mean flow velocity |V| ~ 5.57ms -1 . This result is in good 
concordance with alternative measurements of the velocity (hot wire anemom- 
etry). 

The pertinence of both the tf-rid approach, and the length scale selection, within 
the context of these experiments is illustrated by the high similarity between the 
tf-rid obtained from the two different (and independent) measurement channels. 
Both channels evidence the existence of time localized structures (energy pack- 
ets in the tf-rid) at the same set of instants and with identical Doppler shifts, 
which was to be expected as the probed wave vectors are very close to each 
other (in modulus and direction) and are defined within the same measurement 
volume (cf figure (|l|)). 

The interest of such tf-rid representation in the present context is thus threefold. 
First, it emphasizes the fact that the broadening of the averaged spectrum is 
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due to both the short duration of the structure or events that are the origin 
of the scattering effects, and the fluctuations of the Doppler shifts (the center 
of mass of the energy packets in the tf-rid are not well aligned along an axis 
A.v = constant. Second, such representation provides a straightforward way to 
detect the time of occurrence of events at a given scale. A systematic study 
of the statistical distributions of such events, as functions of the probed scale 
is presently a matter of great interest and will be presented in a future paper. 
Third, the use of two simultaneous channel measurements allows one to test the 
existence of events at different wavelengths arbitrarily close to each other (scale 
resolution). In the preceding example, the characteristics were set to probe two 
close scales (namely |g sca t,i| = 369, 6to -1 and \q SC at,2\ = 371,6m -1 ). Structures 
that are evidenced on one channel clearly appear jointly in the second channel; 
this latter point is illustrated on the central plot of figure (^), where the nor- 
malized time marginal of the geometric mean of the rid of the channels, d(t) is 
shown: 

/ f "r„ 5 ' tfr ridA (t, f n ).tfr rid , 2 (t, f n )df n 

d (t) = I == rT - =5 

V lfn=0 t f r nd,l{t, f n )y J f " =0 tfr rldy2 (t, f n ) 

tfr r id,i is the tf-rid of the signal recorded on channel i, and is the normalized 
frequency. Note that < d(t) < 1 by virtue of Cauchy-Schwarz inequality. This 
method (maximum likelihood criterion) allows detection of the joint presence of 
energy at a given time an d with a given Doppler shift on both channels. Figure 
allows to evidence the similarity of the signals recorded on the channels, 
even when considered in details, thus revealing similar detection of the local 
dynamics. 

The same experiment was performed by resetting the analyzing frequency 
of the second channel to v$ : 2 — 35kHz. Thus, \q SC at,2\ — 442, 4m -1 with the 
same geometry, i.e. q SC at.i and q sca t,2 were still colinear, aligned with the flow 
and oriented downstream. The estimated tf-rids obtained for this situation are 
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shown on figure (g). The normalized co-occurrence function d(t) takes lower 
values than those obtained in the preceding case. Furthermore, the number 
of co-occurrences that are detected is much lower in this case, compared to 
the co-occurrence detected when the tuning condition (cf eq. |8|) is verified. 
Further insight is gained by repeating the previous experiments with \q sca t i| 
held fixed, while |<jf aC ot,2| is varied according to vq 2 ranging from 25kHz to 
35kHz (i.e. 316m" 1 < \q sca t,2\ < 442m -1 ). For each couple (q soa t,i, Qscatfi), 
defining a spatial wavector separation Sq = \q sca t,2 — Qscat,i \ , we compute the RID 
representations. Then we introduce a time-scale generalized cross-correlation 
function C(r, Sq) defined as: 



where (») t stands for an average over the time variable t. Figure (|7|.a) sketched 
the results of this series of experiments on a 2D-plot where the amplitude of 
C(t, 5q) is coded in gray levels, as a function of 5q, and the delay time variable 
r. The evolution of C(r,Sq) for 3 values of Sq (namely Sq — 0.0,0.2,-0.2) is 



time r evolves monotonously from negative values when Sq < towards positive 
values when Sq > 0. The evolution of the maximum of the correlation function is 
plotted on figure (@.d) as a function of Sq. As expected, the correlation function 
is maximum for Sq = decreases smoothly as \Sq\ increases. The observed 
evolutions of both the value and the sign of the delay time is an experimental 
evidence of a time-continuous turbulent cascade of energy from large scales 
(Sq < 0) towards smaller scales (Sq > 0). Up to our knowledge, this time 
continuity of the cascade which is a direct evidence of the time reversal symmetry 
breaking had not been yet experimentally demonstrated. From Figure (0.b) we 
measure a typical extension of the detected wave-packets in the q-domain : 
Sq ~ 0.4cm -1 , corresponding to a spatial correlation length of about 10 cm. 
It is worth noticing here that the measured correlation length remains smaller 




plotted on Figure (0.b). Figure (|?].b and .c) demonstrates that the mean delay 
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than the size of our transducers. It is interesting at this point to mention a 
few words about the spatial and spectral resolutions of the experimental setup. 
The spatial resolution is determined by the maximum frequency one is able to 
propagate through the flow of interest; actually, this frequency is around 200kHz 
enabling to probe lengthscales as small as 1.5mm [^6| (this upper frequency limit 
is due to sound absorption phenomena). The spectral resolution, determining 
the minimum spatial wavevector separation one can probe, depends directly on 
the size L of the transducers which fixed the width Aq of the equivalent spatial 
filter |Q (Aq ~ j- see also section 2). Notice however, that the latter limitation 
(analogous to the Rayleigh criterion used in optics) can be overcome by resorting 
to a deconvolution scheme once the diffraction pattern of each transducer has 
been modelized or characterized. 

6 Conclusion 

We propose a method for evidencing coherent vorticity structures in turbulent 
flows. Though this existence seems to be widely accepted and evidenced by 
numerical simulations Q or in some peculiar flows |3|, || it still deserves further 
experimental efforts. Furthermore, no clear definition seems to exist for co- 
herent structures. The experimental technique described here is derived in the 
spatial Fourier domain and recovers the intuitive notion of phase stationarity 
thus allowing one to understand these structures as propagating energy (more 
precisely enstrophy) packets. 

The main advantage of the acoustic technique relies on its global and non pertu- 
bative character, as it gives access to the instantaneous value of spatial Fourier 
components of the vorticity field. 

Using two distinct channels (acoustic interferometry) we are able to establish 
the robustness (noise immunity) of our vorticity detection scheme. In addition, 
this configuration enables the simultaneous probing of two spatial wave vec- 
tors arbitrarily close to each other (scale resolved spectroscopy). This latter 
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remark is in contrast with usual one-point measurement (hot wire anemometry) 
in turbulence which do not allow a clear separation of time and space dynamics. 
Indeed in such measurements, one has to use the Taylor hypothesis in order to 
recast time increments into spatial increments. 

The joint resort to acoustic scattering and to time-frequency distributions leads 
to a measurement process whereby time (date of occurrence) and frequency 
(Doppler shifts) informations are preserved (time resolved spectroscopy). It is 
worth noticing here that this kind of analysis could be easily transposed to data 
obtained from numerical simulations. Indeed, in the case of pseudo-spectral 
direct numerical simulations (DNS), a great part of the computing time is spent 
in the spatial Fourier space. We claim that time-frequency distributions offer an 
efficient tool for the detection of coherent structures as it is a "blind" detection 
: no a priori information is needed (for example on the shape of the structures) 
and one only resorts to a stationary phase criterion (stating that components 
of a coherent structure are convected with the same velocity). From the RID 
representations of the scattered signals (simultaneously acquired), we define 
a generalized time-scale correlation function based upon a maximum likehood 
criterion. This time-scale correlation function reveals salient features of the 
time-scale dynamics of the detected coherent structures. The most interesting 
one is the existence of a cascade process from larger scale towards smaller one, 
this process being continuous in time. This latter observation indicates that 
coherent vorticity structures play a crucial role in the time reversal symmetry 
breaking, a main feature of the transition to turbulence. 

Finaly, we wish to underline here, the similarities between our experimental re- 
sults and analytical models of coherent structures in the form of spiral vortices 
recently introduced by T.S. Lundgren p3|. We would also like to underline the 



similarity between our experimental approach, and recent numerical 42 and 
theoretical [^lj approaches relying on wavelet decompositions of the turbulent 
velocity field (one-point measurement). 
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FIGURE CAPTIONS 



Figure [l]: Experimental setup. 

Figure ^: Averaged PDS of direct and differentiated recorded data on a jet- 
flow experiment. The overall length of the data field was 2 18 points . The PDS 
is estimated from 2048 points data segments with 1024 data points overlap. 

Figure ||: (a) Estimated spectra from channels 1 and 2; (b) Estimated spectra 
of differentiated signals recorded on channels 1 and 2; (c) Estimated cross- 
coherence function between channels 1 and 2; (d) phase of the cross-channel 
interspectrum (notice that the frequency span is smaller for this last plot). 

Figure |J 2 channels tf-rid estimation for approximate tuning conditions 
(only positive values are represented). 

Figure ^ 2 channels tf-rid estimation for approximate tuning conditions : 
details. 

Figure]^: 2 channels tf-rid estimation for q sca t.i = 369,6m -1 and q sca t,2 = 
442, 4m -1 (only positive values are represented). 

Figure ^[ Evolution of the time-scale correlation function C(r,Sq). (a) 2D 
plot versus time delay r and spatial wavectors separation Sq; (b) Representations 
of C(r,Sq) for three different values of 8q (0.0cm -1 , 0.2cm -1 , —0.2cm -1 ); (c) 
Plot of the mean delay time versus 5q; (d) Plot of the maximum of C (r, Sq) 
versus Sq. 
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31 




Figure 2: 
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Figure 3: 
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